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A photovoltaic (PV) emulator (PVE) is essential equipment for the research 
and diagnostic of PV generation. It is a convenient, highly efficient, and 
low-cost approach when compared to controllable light sources. 
Nonetheless, the implementation of the partial shading capability in a PVE is 
highly limited in terms of efficiency, computation burned, number of power 
converters, and flexibility to change in the ambient condition. This paper 
proposes a partial shading adjuster for a PVE that can overcome the 
aforementioned limitations. The adjuster is applicable to the conventional 
PVE since it is based on an algorithm that can be added to the controller of 
the PVE. By adding the adjuster, the conventional PVE can emulate partial 
shading. The partial shading adjuster is added into a PVE that uses the direct 


PI controller referencing control strategy with the buck controller regulated by the 
PV proportional-integral controller. The results show that the PVE maintains its 


Single diode model accuracy and produces a stable output voltage and current during the load 
changes when the adjuster is added. In conclusion, the proposed partial 
shading adjuster able to improve the capability of the PVE by providing a 
real-time partial shading capability. 
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1. INTRODUCTION 

Photovoltaic (PV) emulator (PVE) has become essential equipment for the research and 
development of the PV generation-based system. The common control strategy for the PVE is the direct 
referencing method. The control strategy consists of the standard PV model (PVM), proportional-integral 
(PI) controller, and buck converter. It is a simple control strategy that does not require any additional 
algorithm to operate since the operating point is determined by the transient response of the PVE. 
Nevertheless, this PVE does not consider the partial shading effect. 

Partial shading occurs when the PV modules do not receive similar irradiance. This is due to 
shadows cast by clouds or buildings or the PV modules is not cleaned properly which leads to the dust layer. 
Computing partial shading is a difficult task especially if it is implemented into the PVE. As a result, the 
majority of the research does not include partial shading capability. There are 2 approaches that can be used 
to emulate partial shading for the PVE, which are the hardware and software approaches. For the hardware 
approach, the amplifier circuit combined with the PV cell is a simple PVE with partial shading capability [1]. 
This method amplifies the current-voltage (I-V) characteristic of the PV cell, which allows high-power 
emulation at a lower cost. Nonetheless, the method is highly inefficient since it operates in the linear region 
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and is limited to the type of PV cell used. Another method is to build the PVE based on the transistor [2]. 
This method uses the PVM electrical circuit model as a reference and the diode presented in the PVM is 
replaced with the transistor. Although this method is simple and it requires a low cost to build, it has several 
limitations. This method requires multiple power sources, has low efficiency, and has a high ripple when 
emulating the certain region on the I-V characterised curve. Another method used to emulate partial shading 
is to combine several PVE [3]. This approach is suitable for a simple series configuration since the digital 
controller used needs to compute multiple PVM in real-time, which burden the digital controller and may 
affect the performance of the PVE. This method is also unable to emulate properly at a certain region of the 
I-V characteristic curve. 

For the software approach, the control strategy is based on the common direct referencing method. 
The modification is done on the PVM to allow partial shading capability. If the PVM output is added 
together based on the PVM input, these PVMs are able to produce partial shading [4], [5]. This control 
strategy is called the PVM adder method. It is nearly similar to the multiple PVE method with the advantage 
of using only a single power converter. However, the digital controller stills need to compute multiple PVM 
to emulate the partial shading. To reduce the computation burden, the computation of the PVMs is conducted 
offline and the data is collected and placed in the look-up table (LUT) for the real-time operation [6]—[9]. 
Although this method solves the computation problem, it has a flexibility problem. This means that it is 
difficult to change the irradiance and temperature during real-time operation since the I-V curve is generated 
externally during the offline period. 

The literature shows that emulating partial shading is a difficult task. When it comes to PVE with 
the partial shading capability, Recently, the single-based PVM with partial shading capability is developed, 
which significantly reduce the processing burden [10]. Therefore, this model is suitable to solve the 
limitation faced by the PVM adder method. This paper presents a new method for the PVE with the ability to 
produce partial shading based on the software approach. It uses a new adjuster that can produce partial 
shading by only adjusting the input of the standard PVM. This PVM is based on the single diode model with 
a series resistor or IDIR model. The common direct referencing method is chosen as the control strategy. 
The continuous current mode buck converter is selected for the PVE and controlled using the common PI 
controller. The PVE is simulated using the MATLAB/Simulink simulation platform. The next section 
discusses the control strategy used in the PVE. Then, the design of the MPPT converter is shown in the 
following section. The results and discussion are covered in the next section. The paper ends with a 
conclusion. 


2. PROPOSED PHOTOVOLTAIC EMULATOR 

The PVM used for the PVE is a simplified version of the PVM adder method. Instead of using 
multiple standard PVM, the irradiance (G), module temperature (T), and PV voltage (Vp) are adjusted so that 
only one standard PVM able to be represented as multiple standard PVM for the partial shading. The 
adjustment algorithm is called the partial shading adjuster, as shown in Figure 1. During the partial shading, 
there are multiple G and T inputs for the PVM (Gec and Tec). The Vp input for the PVM is higher than the 
open-circuit voltage, Voc, of a single PV module since it is Vp in series (Vpy sr). The function of the partial 
shading adjuster is to select a suitable G and T (G; and 7;) for the one standard PVM. Then, the V,, is reduced 
to for the use of the one PVM (Vp, needs to be equal to or lower than Voc for a valid PVM). The new V,, is 
called adjusted Vp, Vpy_aaj. The final output for the PVM is the series PV current, Jp, sr. 

The Ipv sr becomes the reference current, Les, for the PI controller. The Lef is compared with the 
output current, Jv, from the buck converter to obtain the error for the PI controller. The PI controller uses the 
error to produce the corresponding duty cycle, d. The pulse width modulation, PWM, use the d and produce 
the switching pulse, sp, for the buck converter. The buck converter produces the J, according to the Zes. The 
corresponding output voltage, Vo, from the buck converter is sent to the partial shading adjuster to become 
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Figure 1. The block diagram of the proposed PVE with partial shading capability by adding the partial 

shading adjuster 
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2.1. Photovoltaic model 

The PVM consists of two components, which are the standard PVM and the partial shading adjuster. 
The purpose of the standard PVM is to compute a standard I-V characteristic curve without partial shading. 
While the function of the partial shading adjuster is to receive multiple irradiance and temperature input, 
select one irradiance and one temperature, and sent it to a single standard PVM. The selection allows the 
production of the partial shading characteristic. 


2.1.1. Standard photovoltaic model 

The mathematical characteristic equation for the standard PVM is presented using (1), which 
calculate the photovoltaic current, /,, [11], [12]. It is based on the single diode model with a series resistor, 
which is widely used in the PVE application due to the low computation burden [13]—[15]. Commonly, the 
parallel resistor, Rp, is considered in the single diode PVM. Nonetheless, it is not applicable for the partial 
shading adjuster [10]. The Amerisco 80J-B PV module is chosen for the simulation, in which the number of 
PV cells in a module, ns, is 72 [16]. The K» and temperature coefficient short circuit current, Ku, 
are -0.16 V/°C and 0.002 A/°C, respectively. The PV module has the maximum power point voltage and 
current at the STC (Vmp_st and Imp_sic) at 35.8 V and 2.23 A, respectively. The open-circuit voltage and short 
circuit current at the STC (Voc_stc and Isc_stc) are 44.4 V and 2.32 V, respectively. 


Te = Ipn — 1, [exp (ee) -1] (1) 


Where Jpn is the photo-generation current, Vpy is the PV voltage, J; is the saturated current, V, is the thermal 
voltage, and A is the ideality factor. 


2.1.2. Partial shading adjuster 

The partial shading adjuster is an essential part of the PVE to produce partial shading. The algorithm 
starts by loading the parameters of the PV module and constants, as presented in Figure 2. The Jp} for the first 
G and T (G; and Ty) is calculated. Then, the V, and J, are evaluated, respectively. Then, the Jpn for the 
following G and T are calculated until the ns is achieved. At the same time, the critical PV voltage, Vpy cri, iS 
determined using (2) and (3). Then, the sum of PV voltage, Vp sun, is calculated using (4). The Vpy cri is 
compared with the V,, sr. If the Vp sr is smaller than Vpy sum, the Vpy sum is updated by subtracted with Vpy cri 
and the Vy aaj is calculated using (5). This condition occurs when the operating point of the PVE is located 
on the I-V characteristic curve for the highest G. 


Vov_criO = ABV, (T(t) - pn (GCI), T(i)) Rs (2) 
paii p a aD P 1] 3) 
Is(T(i)) 
Vpv_sum = Eia [Vv erii] + (i a 1) ‘aaj Vyv_cri_adj (4) 
Vpv_aaj = Vov_sr = Vov_sum (5) 
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Figure 2. The flowchart of the PVM with the partial shading capability based on the partial shading adjuster 
method 
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If the Vp sr is equal to or larger than Vp cri, the operating point of the PVE is located on the I-V 
characteristic curve for the second-highest and lower G. In this condition, the algorithm checks whether the 
current partial shading adjusted, i, is equal to or lower than the N. If the condition is not true, this means that 
the computation is already at the lowest G and the Vy aaj is calculated using (5). If the condition is true, this 
means that the computation is between the highest G and the lowest G. For this condition, the computation is 
repeated until the Vpv_sr is smaller than Vpy_cri. 

After the computation of the partial shading adjuster is completed, the Vp aq, Gi, and T; are 
produced and used by the standard PVM. Since the standard PVM is based on an implicit equation, the 
iteration method such as the Newton-Raphson method is needed. Using the Newton-Raphson method, (1) is 
iterated until the number of iterations for the standard PVM, J, is achieved. In the end, the J,, becomes the 


Ips Te 


2.2. Buck converter 

The PVE requires a power converter to work. The power converter used for the PVE is commonly 
the buck converter due to the low passive components requirement and high-power efficiency [6], [17], [18]. 
The equivalent circuit of the buck converter provided in Figure 3 contains parasitic resistors for a more 
realistic simulation. There are two parts in designing the buck converter, which is the design of the passive 
element and the controller for the buck converter. 


Figure 3. The topology of the buck converter with parasitic resistors 


The design of the buck converter for the PVE is based on the conventional buck converter [19], [20]. 
The first parameter calculated is the input voltage, V; using (6). The V; is affected by the maximum Voc, 
Voc_max, and maximum d, dmax. Based on the simulation of the PVM, the Voc max is 195.6 V. While the dmax is 
set to 0.9. The calculated V; is 217.3 V. Since the calculation does not consider the nonideality factor, a 
higher voltage needs to be chosen, which is 250 V. 


Voc max 
V, = ema © 


dmax 


The inductance, L, design is important to maintain continuous current operation. This is essential 
especially for the PI controller design since a different controller is needed if the buck converter operates in 
the discontinuous current mode. The Vo max and the maximum output resistance, Ro max, balance needs to be 
considered in the design. In the ideal condition, the Ro maxis infinity. Nonetheless, the buck converter cannot 
be controlled when the Ro maxis infinity [19]. Therefore, a certain limit is chosen, which is 700 Q. This Ro max 
is chosen because it covers the most area in the I-V characteristic curve. Based on the Ro max, the output 
voltage corresponding to Ro max, Vo_Ro_max, 18 obtained from the PVM. Using (7), the required L is 3 mH, 
which the switching frequency, fs, is 20 kHz. The corresponding internal resistance of the inductor, rz, is 
0.9 ©. 


Vo Ro max 
(1-2 Ro max V; Ro max 
i 


L= (7) 


2fs 


The capacitance, C, is designed to operate the buck converter at the desired ripple factor. The 
chosen output voltage ripple factor, Yvo, for the PVE is 0.5% [21]. Based on (8), the C required is 36 uF, 
which minimum D, dnin, is 0.01. The corresponding internal resistance of the capacitor, rc, is 0.3 Q. 


1-dmin 
= 8 
BLYV fs (8) 
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To design the PI controller, the transfer function of the power converter is needed. The transfer 
function of the buck converter, G}(s), is needed, as provided in (9) [18], [22]. The G,(s) has the J, as of the 
output since the current-controlled direct referencing method is used for this PVE. For the transfer function 
of the PI controller, G,(s), it can be presented as (10). Using the automatic tuning function provided by the 
MATLAB/Simulink, the proportional and integral gains (K, and K;) are 0.0001 and 39, respectively. 


1 
= ig(s) a Vi LC 
G(s) = A) Ro saent (9) 
_ âs) _ K; 
Gels) = 3 = Kp + Is (10) 


3. RESULTS AND DISCUSSION 

There is five performance analysis conducted for the proposed PVE. This performance analysis is 
conducted to ensure the direct referencing method control strategy able to cope with the partial shading 
adjuster modification. The accuracy, load, irradiance, stability, and efficiency tests are conducted using a 
resistive load. While the compatibility test is conducted using the MPPT converter as the load for the PVE, 
which the resistive load is removed from the simulation. 


3.1. Accuracy test 

The accuracy is an essential aspect when it comes to PVE. The PVE is considered accurate if the 
operating points fall on the I-V characteristic curve of the PV. Three types of accuracy tests can be conducted 
on the PVE. These can be either by comparing the V, and J, of the PVE with the V, and Zp of the PVM [23], 
[24], PV module [25], or both PVM and PV module [11]. The comparison with the PV module evaluates the 
accuracy of the control strategy and PVM used in the PVE. While the comparison with the PVM only 
evaluates the accuracy of the control strategy used in the PVE. Since the paper focused on the accuracy of the 
control strategy and not the accuracy of the PVM, the V, and J, of the PVE are compared with the V,, and Jp, 
of the PVM. 

The comparison of the I-V characteristic curves produced by the PVE and PVM are shown in Figure 4. 
The result shows that the V, and J, produced by the PVE is on the I-V characteristic curve of PV. This show 
that the proposed PVE has high accuracy. The I-V characteristic of the PVE is in the discrete form because 
the accuracy of the PVE becomes lower if the variable load is sweep from zero to infinity. This is due to the 
slow transient response of the PVE. The PVE requires a certain time to achieve a steady-state. Therefore, the 
V, and I, of the PVE are recorded after the steady-state is achieved with multiple fixed loads. There is also no 
short-circuit (Ro is zero) and open-circuit (Ro is infinity) tests conducted on the PVE. This is because the buck 
converter cannot work properly during this condition. When the R, is equal to zero, the V, needs to be zero, 
which then produce zero J, if a buck converter is used. However, the py from the PVM is not zero when V,, 
equals to zero, which means that the PVE fails to emulate the short-circuit test. If the Ro is infinity, the Vo 
cannot be regulated and operates near to V;. A different type of converter is needed in order to operate at 
these points. 
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Figure 4. The comparison of the I-V characteristic curve of the PVM and PVE 
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3.2. Load test 

The load test is conducted on the PVE to analyse the transience performance and stability during 
load changes. A good PVE can change to a new operating point without any oscillation at a fast rate. By 
referring to Figure 5 (a), the load starts at 30 Q. At 0.2 s, the load is stepped up to 160 Q for 0.2 s. Then the 
load is stepped down back to 30 Q. The waveforms of the Vo and Io of the PVE are recorded in Figure 5 (a) 
and Figure 5 (b), respectively. The percentage overshoot, %OS, of the is observed and it is calculated using 
(11) [26]. A PVE needs to have zero or low %OS to avoid damaging the load or the components in the PVE. 
Another factor that needs to be considered is the settling time, ts, which is used to analyse the transient 
performance of the PVE, which is defined as the time taken for the Io to be within 2% of its steady-state 
value. 


%OS = roma ios x 100% (11) 


o(ss) 


Where Lomax) is the maximum output current overshoot and Iss) is the steady-state output current. 

Based on Figure 5 (b), when the load is stepped down from 160 Q to 30 Q, there is a large I, 
overshoot occur. Using (11), the %OS is 147.7%. The large J, overshoot is caused by the C that try to 
maintain V, after the load is stepped down. This overshoot cannot be avoided when the buck converter is 
used for the PVE. Nonetheless, the C is only 36 uF and the J, should be reduced quickly. However, the Zo 
waveform shows that it takes 43.4 ms to achieve a steady-state. The slow response is due to the constraint of 
the PI controller. To reduce the duration of the overshoot, a faster and more robust controller like fuzzy logic 
or shift controllers [20]. The result shows the ¢, during the beginning of the operation is 44.1 ms. When the 
load is stepped up and down, the f; is 80.5 ms and 43.4 ms, respectively. The results show that the ts is faster 
when the load is stepped down compared to when the load is stepped up. The slow response during a high 
load condition is due to the limitation of the PI controller. The PI controller is designed at the lowest load to 
avoid unstable output. If the PI controller is designed at the higher load condition, the output becomes 
unstable when the load is low. Nonetheless, the PI controller alone does not affect the t, for the PVE [20]. If 
the effect of PI controller alone is accounted, the ts increases as the V, increases. By referring to Figure 6, this 
relationship does not apply. This is because the control strategy used in the PVE affects the t, for the PVE. 
Since the direct referencing control strategy is used for the PVE, the t; becomes higher at the local and global 
MPP. The ft; becomes lower when the operation leaves the MPP and move to the constant voltage region. A 
robust PVE should maintain a constant t, at the various operating point. 
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Figure 5. The time-domain response of the PVE when the load is stepped up from 30 Q to 160 Q at 0.2 s and 
stepped down from 160 Q to 30 Q at 0.4 s (a) output voltage waveform (b) output current waveform 
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Figure 6. The settling time of the output voltage of the PVE at various load condition 
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4. CONCLUSION 

The PVM requires a high computation since it is based on an implicit equation that requires multiple 
iterations. To produce partial shading, several PVMs need to be used, which further increases the 
computation burden. Nonetheless, the partial shading adjuster allows partial shading capability with only a 
single PVM. This gives an advantage for the PVE because the PVE able to operate in real time and requires a 
low computation to work properly. The partial shading adjuster is a new concept and is not being tested in the 
PVE. Based on the simulation results, this adjuster is suitable for the PVE. The accuracy test shows the PVE 
can imitate the I-V characteristic curve produced by the PVM. This result not only shows combining the 
partial shading adjuster with the PVE has produced high accuracy results but also shows it is stable during 
the steady state operation. The result shows that the transient response becomes slower as the load increases. 
There is also a current spike that occurs during the load stepped down. A different power converter, closed- 
loop controller, and control strategy are needed to improve the PVE. In conclusion, the partial shading 
adjuster can work with the PVE since it is accurate. The transient response of the PVE requires improvement 
due to the slow response at a high load and a high current spike when the load decreases. 
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